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I. INTRODUCTION 

The idea that spacetime might emerge from more fundamental degrees of freedom has 
long fascinated physicists. The holographic principle suggests that a (H + 1)-dimensional 
spacetime might emerge from degrees of freedom in a H-dimensional theory without gravity 
[1, 2]. While a completely general implementation of this idea is still lacking, the AdS/CFT 
correspondence provides a specihc example in which to probe the holographic emergence 
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of spacetime. AdS/CFT is a conjectured correspondence between D-dimensional conformal 
field theories (CFTs) in Minkowski space and {D + l)-dimensional asymptotically anti-de 
Sitter (AdS) spacetimes [3-5]. An intriguing aspect of this duality is the Ryu-Takayanagi 
formula [6, 7], according to which the entanglement entropy of a region B on the boundary 
is proportional to the area of a codimension-two extremal surface B embedded in the bulk 
curved spacetime whose boundary is B: 

aTea(B) 

S{B) = — ^ ^ -I- corrections. (1) 

^ ' AG 

In other words, given a CFT state, one may think of bulk distance and geometry (at least 
near the boundary) as being charted out by the entanglement properties of the CFT state. 

A central question in this picture of spacetime emerging from entanglement is: What is 
the precise relationship between bulk degrees of freedom and boundary degrees of freedom? 
Expressed in a different way, what is the full map between states and operators in the 
boundary Hilbert space and those in the bulk? While investigations of AdS/CFT have 
thrown a great deal of light on this question, explicit simple models are still very helpful for 
studying it in more detail. 

Meanwhile, from a very different perspective, tensor networks have arisen as a useful way 
to calculate quantum states in strongly-interacting many-body systems [8]. One signihcant 
example is the Multi-scale Entanglement Renormalization Ansatz (MERA) [9], which is 
relevant for critical (gapless) systems, i.e., CFTs. Starting from a simple state in a low¬ 
dimensional Hilbert space, acting repeatedly with hxed tensors living on a network lattice 
produces an entangled wave function for the quantum system of interest; varying with respect 
to the tensor parameters efficiently computes the system’s ground state. 

Working “backwards” in the MERA, starting with the ground state and gradually re¬ 
moving entanglement, produces a set of consecutively renormalized quantum states. This 
process reveals a renormalization direction along the graph, which may be thought of as an 
emergent radial direction of space. As pointed out by Swingle [10], the MERA graph can 
serve as a lattice discretization of spatial slices of AdS. Furthermore, one can use the MERA 
to calculate the entanglement entropy of regions of the original (boundary) critical system; 
this calculation amounts to tracing over bonds in the tensor network that cross the causal 
cone of the boundary region. The causal cone is a sort of extremal surface for the MERA, 
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motivating comparison to the Ryu-Takayanagi formula. 

It is therefore natural to conjecture that the MERA provides a concrete implementation 
of the emergence of spacetime, in the form of a correspondence between boundary and 
bulk regions reminiscent of AdS/CFT [10]. Such an AdS/MERA correspondence would be 
extremely useful, since the basic building blocks of the MERA are discrete quantum degrees 
of freedom from which quantities of physical interest may be directly calculated. Some 
specihc ideas along these lines have recently been investigated [11-14]. 

In this paper, we take a step back and investigate what it would mean for such a cor¬ 
respondence to exist and the constraints it must satisfy in order to recover properties we 
expect of physics in a bulk emergent spacetime. After reviewing the MERA itself and pos¬ 
sible construals of the AdS/MERA correspondence in the next section, in Sec. Ill we then 
derive relationships between the MERA lattice and the geometry of AdS. We find that the 
MERA is unable to describe physics on scales shorter than the AdS radius. In Sec. IV we 
explore constraints from calculating the entanglement entropy of regions on the boundary, 
in which we are able to relate MERA parameters to the central charge of the GET. Finally, 
in Sec. V we apply the covariant entropy (Bousso) bound to regions of the bulk lattice. In 
the most naive version of the AdS/MERA correspondence, we hnd that no combination of 
parameters is consistent with this bound, but we suggest that generalizations of the tensor 
network may be able to provide a useful correspondence. 

II. ADS/MERA 

Let us begin by recalling the dehnition and construction of the MERA. We will then 
introduce the AdS/MERA correspondence and discuss the motivation for and consequences 
of this proposal. 


A. Review of the MERA 

The MERA is a particular type of tensor network that provides a computationally efficient 
way of Ending the ground states of critical quantum many-body systems, i.e. CFTs, in D 
dimensions. (For a recent review of tensor networks in general, see Ref. [8]. Detailed analyses 
of the MERA are given in [9, 15, 16] and references therein.) In this work, we restrict our 
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FIG. 1. (a) Basic construction of a /c = 2 MERA (2 sites renormalized to 1). (b) The squares 
represent disentanglers: unitary maps that, from the moving-upward perspective, remove entan¬ 
glement between two adjacent sites, (c) The triangles represent isometries: linear maps that, again 
from the moving-upward perspective, coarse-grain two sites into one. Moving downward, we may 
think of isometries as unitary operators that, in the MERA, map a state in E (8) |0) into V 0V. 
The i and j labels in (b) and (c) represent the tensor indices of the disentangler and isometry. 


attention to the case D = 1 -|- 1. 

The MERA tensor network is shown in Fig. 1. The qnantnm system being modeled by 
the MERA lives at the bottom of the diagram, henceforth “the bonndary” in anticipation of 
the AdS/MERA connection to be explored later. We can think of the tensor network as a 
qnantnm circnit that either rnns from the top down, starting with a simple inpnt state and 
constrncting the bonndary state, or from the bottom np, renormalizing a bonndary state via 
coarse-graining. One dehning parameter of the MERA is the rescaling factor fc, dehning the 
nnmber of sites in a block to be coarse-grained; in Fig. 1 we have portrayed the case k = 2. 
The sqnares and triangles are the tensors: mnltilinear maps between direct prodncts of vector 
spaces. Each line represents an index i of the corresponding tensor, ranging over values from 
1 to the “bond dimension” y. The boundary Hilbert space "Hboundary = E®-^boundary jg given 
by a tensor product of A/boundary individual spaces V, each of dimension y. (In principle 
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the dimension of the factors in the boundary could be different from the bond dimension of 
the MERA, and indeed the bond dimensions could vary over the different tensors. We will 
assume these are all equal.) 

As its name promises, the MERA serves to renormalize the initial boundary state via 
coarse-graining. If we were to implement the MERA for only a few levels, we would end 
up with a quantum state in a smaller Hilbert space (dehned on a hxed level of the tensor 
network), retaining some features of the original state but with some of the entanglement 
removed. However, we can also run the MERA backwards, to obtain a boundary state from 
a simple initial input. By varying the parameters in the individual tensors, we can look for 
an approximation of the ground state of the GET on the boundary. Numerical evidence 
indicates that this process provides a computationally efficient method of constructing such 
ground states [16, 17]. 

The tensors, or gates, of the MERA come in two types. The first type are the disen- 
tanglers, represented by squares in Fig. 1. These are unitary maps U: V®V^V®V, as 
in Fig. lb. The name comes from thinking of moving upward through the network, in the 
direction of coarse-graining, where the disentanglers serve to remove local entanglement; as 
we move downward, of course, they take product states and entangle them. The second 
type of tensors are the isometries, represented by triangles. From the moving-downward 
perspective these are linear maps W \ V ^ V ®V] moving upward, they implement the 
coarse-graining, see Fig. Ic. The isometries are subject to the further requirement that 
W^W = Jy, where Jy is the identity map on V, and WW^ = Fa, where Pa is a projector 
onto some subspace A <Z V®V. From the top-down perspective, we can also think of the 
isometries as bijective unitary operators Wjj : V ®V ^ V®V, for which a fixed “ancilla” 
state (typically the ground state |0)) is inserted in one of the input factors, as shown in 
Fig. Ic. More generally, isometries could map q < k sites onto k sites, W : —)■ . 

The MERA is not the simplest tensor network which implements coarse-graining. For 
instance, the tree tensor network [18] (also considered in a holographic context in Ref. [11]), 
similar to MERA but without any disentanglers, also implements coarse-graining. How¬ 
ever, tensor networks without disentanglers fail to capture the physics of systems without 
exponentially-decaying correlations, and consequently cannot reproduce a GET ground state. 

An example that invites analysis with a MERA is the transverse-field Ising model [19]. 
In 1 -|- 1 dimensions, the model describes a chain of spins with nearest-neighbor interactions 
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subject to a transverse magnetic field. Its Hamiltonian is 


H 







( 2 ) 


where af and af are Pauli operators and where J and h set the strength of the nearest- 
neighbor interactions and the magnetic field, respectively. Notably, the system achieves 
criticality at J = h, where a quantum phase transition occurs between ordered {J > h) 
and disordered (J < h) phases. In this example, the open legs at the bottom of the MERA 
describe the state of the one-dimensional lattice of spins. A single application of disentanglers 
and isometries can be thought of as a true real-space renormalization, producing a lattice 
of spins that is less dense than the preceding lattice by a factor of q/k. 

In general, much information is required to describe an arbitrary MERA. In principle, the 
Hilbert spaces, the disentanglers, and the isometries could all be different. Also, for /c > 2, 
there is no canonical way of laying out the disentanglers and isometries; the circuit itself must 
be specified. We will restrict ourselves to the case g = 1, so that isometries have 1 upward- 
going leg and k downward-going legs. Further, without loss of generality, we take the same 
vector spaces, disentanglers, and isometries everywhere in the MERA, a simplification that 
is enforced by the symmetries of the boundary ground state. These symmetries — namely, 
translation- and scale-invariance — dictate that the MERA parameters and structure be 
homogeneous across the whole tensor network. 

For geometric considerations, it is useful to abstract away all of the information about 
unitary operators and to draw a MERA as a graph as shown in Fig. 2. In such a graph, 
we only indicate the connectivity of sites at any given level of coarse-graining as well as the 
connectivity of sites under renormalization group flow. 

B. An AdS/MERA Correspondence? 

The possibility of a correspondence between AdS and the MERA was first proposed 
by Swingle in Ref. [10], where it was noted that the MERA seems to capture certain key 
geometric features of AdS. At the most basic level, when viewed as a graph with legs of fixed 
length, a MERA may be thought of as a discretization of the hyperbolic plane, which is a 
spatial slice of AdSs. In this discretization, the base of the MERA tree lies on the boundary 
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FIG. 2. (a) A k = 2 MERA, and (b) the same MERA with its disentanglers and isometries 
suppressed. The horizontal lines in the graph on the right indicate lattice connectivity at different 
renormalization depths, and the vertical lines indicate which sites at different depths are related 
via coarse-graining due to the isometries. Each site, represented by a circle, is associated with a 
Hilbert space V with bond dimension y. In the simplest case, a copy of the same Hilbert space 
is located at each site. When assigning a metric to the graph on the right, translation and scale 
invariance dictate that there are only two possible length scales: a horizontal proper length Li and 
a vertical proper length L 2 . 


of the AdS slice and the MERA lattice sites £11 out the bulk of the slice [10, 20]. 

Interestingly, the structure of a MERA is such that it seems to go beyond a simple dis¬ 
cretization of the hyperbolic plane. Certain discrete paths in the MERA naturally reproduce 
geodesics of the hyperbolic plane [10, 21]. Moreover, this phenomenon makes it possible to 
understand the computation of GET entanglement entropy using a MERA as a discrete 
realization of the Ryu-Takayanagi formula [22]. These and other examples [10, 21] seem to 
suggest that a MERA may in fact be elucidating the structural relationship between physics 
on the boundary of AdS and its bulk. 

In this work we take the term “AdS/MERA correspondence” to mean more than simply 
a matching of graph geometry and continuous geometry. In the spirit of the AdS/CFT cor¬ 
respondence, we suppose that (at least some aspects of) both boundary and bulk physics are 
described by appropriate Hilbert spaces "Hboundary and "Hbuik respectively, which must have 
equal dimensions. A full AdS/MERA correspondence would then be a specification of these 
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Hilbert spaces, as well as a prescription which makes use of the MERA to holographically 
map states and operators in ^boundary to corresponding states and operators in Tfbuik and 
vice-versa. To preserve locality in the bulk and the symmetries of AdS, it is natural to 
identify "Hbuik with the tensor product of individual spaces 14uik, each located at one site 
of the MERA. If it exists, this correspondence provides a formulation of bulk calculations 
in terms of the MERA. An AdS/MERA correspondence should allow us to, for example, 
calculate bulk correlation functions, or bulk entanglement entropies using tools from or the 
structure of the MERA. 

There is one straightforward way to construct such a map ^boundary ^ "Rbuik- We have 
noted that the isometries W -.V ^ V®V can be thought of as unitaries Wu \V®V^V®V 
by imagining that a hxed ancillary state |0) is inserted in the hrst factor; for a fc-to-one 
MERA, one would insert k — 1 copies of the |0) ancilla at each site to unitarize the isometries. 
From that perspective, running upwards in the tensor network provides a map from the 
MERA ground state on the boundary to a state | 0 )®(^“bMuik ^ ^®(fc-i)A4iuik^ where at 
each isometry there is a copy of and Abuik denotes the number of bulk lattice 

sites, excluding the boundary layer. As we ultimately show in Sec. V, one has A/boundary = 
{k — l)Abuik- We can then identify "Hboundary = "Hbuik = R^Woundary think of the tensor 
network as a quantum circuit providing a map between arbitrary states Hboundary —t Hbuik- 
In this construction, the MERA ground state on the boundary gets mapped to the factorized 
bulk state | 0 )®(^“^Wbuik^ other boundary states will in general produce entangled states 
in the bulk (keeping the tensors themselves hxed). 

Something very much like this construction was proposed by Qi [11], under the name 
“Exact Holographic Mapping” (EHM). That work examined a tensor network that was not 
quite a MERA, as no disentanglers were included, only isometries. As a result, while there 
is a map Hboundary “t Hbuik, the bouudary state constructed by the tensor network does not 
have the entanglement structure of a GET ground state. In particular, it does not seem 
to reproduce the Ryu-Takayanagi formula in a robust way. Alternatively, we can depart 
from Qi by keeping a true MERA with the disentanglers left in, in which case the bulk 
state constructed by the quantum circuit has no entanglement: it is a completely factorized 
product of the ancilla states. Such a state doesn’t precisely match our expectation for what a 
bulk ground state should look like, since there should be at least some entanglement between 
nearby regions of space. 
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Therefore, while it is relatively simple to imagine constructing a bulk Hilbert space and a 
map between it and the boundary Hilbert space, it is not straightforward to construct such 
a map that has all of the properties we desire. It might very well be possible to find such a 
construction, either by starting with a slightly different boundary state, or by adding some 
additional structure to the MERA. 

For the purposes of this paper we will be noncommittal. That is, we will imagine that 
there is a bulk Hilbert space constructed as the tensor product of smaller spaces at each 
MERA site, and that there exists a map "Hboundary —^ "Hbuik that can be constructed from the 
MERA, but we will not specify precisely what that map might be. We will see that we are 
able to derive bounds simply from the requirements that the hypothetical correspondence 
should allow us to recover the properties we expect of bulk physics, including the background 
AdS geometry and features of semiclassical quantum gravity such as the Bousso bound on 
bulk entropy. 


III. MERA AND GEOMETRY 

If a MERA is a truly geometrical object that describes a slice of AdS, then the graph 
geometry of a MERA should give the same answers to geometric questions as the continuous 
geometry of a slice of AdS. Here, we reconsider the observation by Swingle [10, 21] that 
certain trajectories on the MERA coincide with trajectories in AdS and we investigate the 
constraints that this correspondence places on the graph metric of the MERA. We hnd that 
a MERA necessarily describes geometry on super-AdS length scales, moreover, there is no 
redehnition of the MERA coordinates that results in the proper distance between MERA 
sites mapping to any sub-AdS length scale. 


A. Consistency conditions from matching trajectories 

In order to speak of graph geometry, one must put a metric on the MERA graph, he., 
one must assign a proper length to each bond in the graph of Fig. 2. Presumably, the metric 
should originate from correlations between the sites in the MERA. In the absence of an 
explicit identification of the origin of the graph metric, however, at least in the case of a 
MERA describing the ground state of a GET, it is sensible to identify two length scales. 
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Explicitly, we must assign a proper length Li to horizontal bonds and a proper length L 2 
to vertical bonds. Indeed, translational and conformal invariance guarantee that these are 
the only two length scales in any graph metric one can assign to a MERA for which an 
AdS/MERA correspondence exists. In particular, the ground state of a GET is translation 
invariant, so each horizontal bond in the hnest (UV-most) lattice should have the same 
proper length so as to respect this symmetry. Self-similarity at all scales then requires that 
any horizontal bond at any level of renormalization have this same proper length. There is 
no a priori reason why the vertical bonds should share the proper length of the horizontal 
bonds and indeed we will see that their proper length will be different. However, again 
by self-similarity and translation invariance, all vertical bonds must be assigned the same 
proper length. 

The observation in Ref. [10] that certain paths in the MERA graph coincide with cor¬ 
responding paths in slices of AdS is what established the possibility of an AdS/MERA 
correspondence. Here we will carefully examine these paths and determine what constraints 
the requirements that they match place on MERA parameters, he., on the bond lengths Li 
and L 2 and on the rescaling factor k. 

Consider a constant-time slice of AdSa with the following metric: 


ds2 


(dz^ + dx^). 


(3) 


We will compare the proper lengths of straight horizontal lines and geodesics in the AdS 
slice to the proper lengths of the corresponding paths in the MERA graph. In the AdS slice, 
let 7 i be a straight horizontal line (d^ = 0 ) sitting at z = Zq with coordinate length xq. 
Let 72 be a geodesic whose endpoints lie near the boundary z = 0 and are separated by a 
coordinate distance xq at the boundary. In this choice of coordinates, such a geodesic looks 
like a semicircle (see Fig. 3). It is a straightforward computation to show that the proper 
lengths of these curves are 



IhiUds = —Xq and 
zo 


Note that there is a UV cutoff at z = a <C xq and that we have neglected terms of order 
a/xQ. 
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FIG. 3. A horizontal line ( 71 ) and a geodesic ( 72 ) in a spatial slice of AdSa. 


We fix Li and L 2 by comparing 71 and 72 to horizontal lines and “geodesics” in the MERA, 
respectively. Consider two sites in a horizontal lattice at depth m {i.e., m renormalizations 


of the UV-most lattice) and separated by a coordinate distance Xq in the coordinate system 
shown in Fig. 2. By hat, this lattice sits at Zq = k^a. The number of bonds between the 
two sites at depth m is XQ/{k"^a) (see Fig. 2 for the case k = 2). It follows that the proper 
length of the line connecting the two points is just 


|7i|mera — Li ■ (number of bonds between endpoints) 


(5) 



To have |7i|Ads = |7i|mera, we should therefore set Li = L. 

Similarly, consider two lattice sites on the UV-most lattice separated by a coordinate 
distance Xq. If we assume that Xq ^ a, then the shortest path (geodesic) in the MERA 
connecting the two lattice sites is the path that goes up in the renormalization direction 
and then back down again. The two sites are separated by xo/a bonds on the UV-most 
lattice, so log^(a:o/a) renormalization steps are needed to make the sites either adjacent or 
superimposed. This means that the geodesic that connects the endpoints is made up of 
2\ogf.{xQ/a) bonds (as we have to go up and then back down again, giving the factor of 2 ). 
It follows that the proper length of the geodesic is 


I 72 IMERA — L 2 ■ (number of bonds in the geodesic) 



( 6 ) 


To have | 72 |Ads = | 72 |mera, we should therefore set L 2 = L\n. k. 
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B. Limits on sub-AdS scale physics 


One aspect of the matching of geodesics that is immediately apparent is that the MERA 
scales Li and L 2 that parametrize the proper distance between lattice sites are of order the 
AdS scale L or larger, as was also noted in Refs. [10, 20]. This rnns connter to the typical 
expectation that, in a discretization of spacetime, one expects the grannlarity to be apparent 
on the UV, rather than the IR, scale. That is, snb-AdS scale locality is not manifested in 
the MERA constrnction and must be encoded within each tensor factor [21]. 

One could try to evade this difficulty by attempting to redehne the MERA coordinates 
(x, (those of Fig. 2) as functions of the AdS coordinates (a;, (those of Fig. 3) 

and taking a continuum limit; above, we assumed that the two sets of coordinates were 
simply identihed. That is, suppose = /(x^'^®) and (For example, 

one could consider f{x) = ex for small e and imagine taking the continuum limit, with 
the aim of making Li much smaller than the AdS scale.) If a is still the UV cutoff on the 
AdS side, then in the MERA we have /(a) as the UV-most lattice spacing and g{a) as the 
UV cutoff in the holographic direction. Consider the computation of | 7 i|. From the AdS 
side, we have | 7 i|Ads = On the MERA side, the number of sites spanned by 

^MERA ^ yj'^Ads^ jg f {a), while the holographic coordinate is = k'^g{a). 

Hence, 


1 7l I MERA 


' /(a) g{z^<^^y 


(7) 


Equating lyiUds = |7i|mera = | 7 i|, we have 


^(7 


AdS\ 


d 




I71I = Li 


fix, 


AdS^ 


/(a 




Ads', 


AdS 


( 8 ) 


Since the right side of the hrst equality only depends on x^'^® and the second equality only 
depends on but we can vary both parameters independently, both expressions must 

be independent of both AdS coordinates. Hence, we must have /(x) = e^x and g{z) = e^z 
for some constants and Plugging everything back into Eq. (7) and comparing with 
IbiUds, we again End that Li = L, so no continuum limit is possible. Similarly, in computing 
I 72 I, we note that the number of bonds between the endpoints on the UV-most lattice 
level is x^^^^//(a), so the geodesic connecting the endpoints has 2 log^(x^^^^/e: 3 ;a) bonds. 
On the other hand, we have | 72 |Ads = 2Lln(x^‘^®/a) = 2Lln(x^^^^/£a;a). That is, in 
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equating | 72 |Ads and | 72 |mera, we must again set L 2 = Link. We thus also find that no 
continuum limit is possible in the holographic direction. That is, we have shown that there 
is a constant normalization freedom in the definition of each of the coordinate distances on 
the AdS and MERA sides of any AdS/MERA duality, but such a coordinate ambiguity is 
unphysical and does not allow one to take a continuum limit. One still finds that the physical 
MERA parameters Li and L 2 are AdS scale. This means that there truly is no sense in 
which a discrete MERA can directly describe sub-AdS scale physics without the addition of 
supplemental structure to replace the individual tensors. This fact limits the ability of the 
MERA to be a complete description of the gravity theory without such additional structure. 
It might be the case that one needs a field theoretic generalization of the MERA, such as 
continuous MERA (cMERA) [23-25] or some local expansion of the individual tensors into 
discrete tensor networks with a different graph structure to describe sub-AdS physics, but 
such a significant generalization of the tensor network is beyond the scope of this work and 
in any case would no longer correspond to a MERA proper. 


IV. CONSTRAINTS FROM BOUNDARY ENTANGLEMENT ENTROPY 

Because the MERA can efficiently describe critical systems on a lattice, quantities com¬ 
puted in the MERA on scales much larger than the lattice spacing should agree with GET 
results. In this section, we will compute the entanglement entropy of io contiguous sites in 
the MERA and exploit known GET results to obtain constraints on the properties of the 
MERA. In particular, we will find an inequality relating the MERA rescaling factor k and 
bond dimension y to the GET central charge c. This constraint is interesting in its own 
right, but it will prove critical in the next section when we begin to compute bulk properties. 


A. MERA and GET Entanglement Entropy 

For a (l-l-l)-dimensional GET in a pure state, the von Neumann entropy of a finite interval 
B, which is typically referred to as the entanglement entropy, is known to be [26, 27] 

5(R) = ^ln4, (9) 
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where the length of the interval is much smaller than the system size. Here, is the length of 
the interval in units of the UV cutoff. In the notation of the last section, we have = x^/a. 
In the special case that the CFT is dual to AdS in 2 + 1 dimensions, the central charge is 
set by the Brown-Henneaux formula [28], 


c = 




( 10 ) 


Also note that the length of the geodesic that connects the two ends of B (the curve 72 in 
Fig. 3) is given in Eq. (4) by I 72 I = 2Llnf'o- The Brown-Henneaux relation allows us to 
reproduce the Ryu-Takayanagi formula [ 6 , 29] from the entanglement entropy, 

= ( 11 ) 

where R = 72 is the extremal bulk surface with the same boundary as B. For a boundary 
with one spatial dimension and a bulk with two spatial dimensions, any simply-connected 
region B is an interval, the extremal bulk surface is a geodesic, area(i?) is a length, and G 
has mass dimension — 1 . 


The MERA calculation of the entanglement entropy of sites in the CFT has an anal¬ 
ogous geometric interpretation. Suppose one is given the MERA representation of a lattice 
CFT ground state, he., one uses a MERA to generate the CFT state. Denote by S'MERA(f'o) 
the entanglement entropy of the resulting state restricted to sites. In Ref. [22], it is shown 
that for a specihc, optimal choice of sites, for £0 parametrically large, the following bound 
is placed on S'mera(^o) for a MERA with k = 2: 


5'mera(^o) < 21og2^ofoX- 


( 12 ) 


Parsing the equation above, this bound essentially counts the number of bonds that the 
causal cone of the io sites in question crosses (~ 2 log 2 £ 0 ) and In 7 is the maximum en¬ 
tanglement entropy that a single bond can possess when the rest of the MERA is traced 
out. 

The causal cone of a region B consisting of £0 contiguous UV sites in a MERA resembles 
a bulk extremal surface for the boundary region B. Given £q sites in the UV, their causal 
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cone is defined as the part of the MERA on which the reduced density matrix (or in other 
words, the state) of B depends. An example of a causal cone is illustrated in Fig. 4. 



FIG. 4. Causal cone (shaded) for a set of £o = 6 sites in a MERA with k = 2. The width of the 
causal cone at depth m is ii = 4, (.2 = 3, £3 = 3, f '4 = 3, etc. The crossover scale for this causal 
cone occurs at ?fi = 2. Between the zeroth and first layer, = 2 bonds are cut by the causal 
cone. Similarly, = 2, = 3, etc. 


In particular, note that the number of bonds that a causal cone crosses up to any hxed 
layer scales like the length of the boundary of the causal cone up to that layer. It is in this 
sense that Eq. (12) is a MERA version of Ryu-Takayanagi. Also note that the width of the 
causal cone shrinks by a factor of ~ 1/k after every renormalization step until its width is 
comparable to k. As such, if one denotes the width of the causal cone at a layer m by 
then ^rn is roughly constant for all m greater than some fh (see Fig. 4). The scale fh is called 
the crossover scale. 

For general fc, it is also possible to formulate a bound similar to Eq. (12) for the entan¬ 
glement entropy of £0 sites. For parametrically large we hnd that 


5'mera(4; B) < 4(fc - 1) logfc 4 In X. (13) 

We demonstrate this bound in App. A using techniques that are similar to those developed 
in Ref. [22] . In particular, note that we do not allow ourselves to choose the location of the 
£0 sites in question. As such, we remind ourselves that Amera can depend on the location of 
the region B (and not only its size) by including it in the argument of Smera- This is also 


16 







the reason why our Eq. (13) is more conservative than the optimal bound given in Eq. (12). 


B. Constraining 5mera 

Let us examine Eq. (13) a bit more closely. As discussed in App. A, 4(fc — 1) is an upper 
bound on the number of bonds that the causal cone could cut at any given depth m below the 
crossover scale rfi. (The crossover scale fh is attained after roughly log^ renormalization 
steps.) For a given causal cone, he., for sites at a given location with respect to the 
MERA, let us parametrize our ignorance by writing 


5'mera( 4;^) < 4/ij(fc)logfc41nx, (14) 

where /s(fc) grows no faster than {k — 1) and counts the (average) number of bonds cut by 
the causal cone at any depth up to the crossover scale. Explicitly, 

.. fh—l 

fi’W ^ i- E < ■ (15) 

m=0 

where denotes the number of bonds that the causal cone cuts at the level. 

Each cut bond contributes at most In y to the entropy (the case of maximal entanglement). 
As such, it is instructive to introduce a parameter rjB G [0,1] that describes the degree of 
entanglement of the sites in the causal cone. In doing so we may rewrite the inequality (14) 
as an equality: 

5'mera(4; B) = Afsik) logfc 4 ■ Vb In y. (16) 

The quantity T^slny is the average entanglement entropy per cut bond in the causal cone 
of B. Equivalently, Eq. (16) may be taken as the dehnition of tje- 

This dehnition of of course depends on the location of B and only applies to bonds that 
are cut by the causal cone of B. In what follows, it will be advantageous to have a notion of 
average entanglement entropy per bond that applies to all bonds in the MERA. To this end, 
start with a lattice consisting of Rot sites in total and consider the limit in which the size of a 
region B is unbounded but where the ratio ^o/4ot is held constant (so that B does not grow 
to encompass the whole domain of the GET). In this limit, S'mera(4;-B) —>■ 5'mera(4) and 
/s(fc) —t f{k) should be independent of the exact location of B, i.e., S'mera should exactly 
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agree with Eq. (9). Let us consequently define the average entanglement entropy per bond 
in the MERA: 


r] In y = lim 


5'mera(^o) 


io^oc 4/(A;)logfc(4) ’ 
The quantity rj is then a property of the MERA itself. 


( 17 ) 


Intuitively, one would not expect each individual bond in the MERA to be maximally 
entangled and so it should be possible to constrain r] more tightly than rj < 1. This expec¬ 
tation is made more precise via the following considerations. To begin, consider a MERA 
with k = 2 and examine a pair of isometries at a fixed depth m. As indicated in Fig. 5a, let 
P 2 denote the density matrix of the bonds and ancillae emanating from the two isometries 
and let pi denote the density matrix of the four highlighted bonds below the isometries. We 
again assume that the ancillae are initialized to the pure product state composed of factors 
of |0). Taking into account the ancillae, or in other words promoting the isometries to uni- 
taries, we see that pi and p 2 are related by a unitary transformation, so S{pi) = S{p 2 ). By 
assumption, the state of each ancilla is | 0 ), so p 2 = P 2 < 8 ) | 0 )( 0 | ( 8 ) | 0 )( 0 | for some density matrix 
P 2 . This in turn implies that S{p 2 ) = S{p 2 ) < 21ny. From the definition of p above, the 
entanglement entropy of a single bond is asymptotically given by ^^Iny, so S{pi) ~ 4plny. 
It therefore follows that p < 1/2. 




FIG. 5. A pair of isometries with their ancillae explicitly indicated for a MERA with (a) k = 2 
and (b) general k. The thick bonds below the isometries, the state of which is denoted by pi, 
are unitarily related to the bonds that exit the isometries and the ancillae, the state of which is 
denoted by p 2 - 
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For general fc, the argument is nearly identical. We again begin by considering a pair of 
isometries at a given level m (see Fig. 5b). Analogously with the k = 2 case, let p 2 denote the 
density matrix of the two bonds and 2k — 2 ancillae emanating from the two isometries and 
let Pi denote the density matrix of the 2k highlighted bonds below the isometries. There is 
only one disentangler that straddles both of the isometries in question for any layout of the 
MERA. As such, at most k of the lower bonds enter a disentangler from below and the rest 
directly enter the isometries. Here as well pi and p 2 are related by a unitary transformation 
so that S{pi) = S{p 2 ). Similarly, p 2 = P 2 ® (|0)(0|)®^^“^ for some density matrix p 2 , so 
S{p 2 ) = S{p 2 ) < 21nx. The region described by pi always consists of 2k bonds, so we may 
again asymptotically write *S'(pi) ~ 2kri\ax- It therefore follows that kp < 1, and since 
f{k) < {k — 1), we may write 

rif(k) < . (18) 

We note that, in computational practice, one typically does not use the “worst-case scenario” 
construction explored in App. A; a more conventional construction would result in a tighter 
bound on f{k) and hence a stricter inequality than Eq. (18). For our purposes, however, we 
will remain as conservative as possible and therefore use the inequality (18) in our subsequent 
bounds. 


C. Matching to the CFT 


Finally, we obtain a constraint on fc, y, and p in terms of the central charge c by collecting 
the results of this section. Let us work in the limit where the interval is much larger than 
the lattice spacing, log^^o 3> 1. We have seen that this is precisely the regime in which p 
and f{k) are well-dehned quantities independent of the choice of B. It is also the regime 
in which we can equate the CFT entropy S{io) = (c/3) In^o with the MERA entropy (16). 
Doing so, the central charge is given by 


3L 



12pf{k) 


In k 


(19) 


Then in light of Eq. (18), we hnd that 


c < 12 


k-l \ 
kink J 


Iny. 


( 20 ) 
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To recapitulate, given a CFT with central charge c and a MERA representation of its 
ground state, a necessary condition for a consistent AdS/MERA correspondence is that the 
MERA parameters k and y satisfy the constraint (20). Importantly, this implies that, for 
a well-dehned semiclassical spacetime (for which c S> 1), the bond dimension y must be 
exponentially large in the size of the AdS scale compared to the Planck scale. 

Let us also note that we can still obtain a bound from Eq. (19), albeit a weaker one, with¬ 
out using the result of Eq. (18). Recall that this latter result relies on having unentangled 
ancillae in the MERA. This is not necessarily the case for other tensor network bulk con¬ 
structions, as we will subsequently discuss. As such, if we disregard the result of Eq. (18), 
we still have by virtue of their dehnitions that f{k) < k — 1 and rj < 1. The following 
weaker but more general bound on the central charge therefore follows from Eq. (19) for 
such generalized tensor networks: 

V. CONSTRAINTS FROM BULK ENTANGLEMENT ENTROPY 

In addition to the compatibility conditions from geodesic matching and boundary en¬ 
tanglement entropy, it is well-motivated to seek out any other possible quantities that can 
be computed in both the MERA and AdS/CFT frameworks, so as to place further con¬ 
straints on any AdS/MERA correspondence. One important example of such a quantity is 
the entropy associated with regions in the bulk, as opposed to on the boundary. 


A. The Bousso Bound 

The notion of placing bounds on the entropy of regions of spacetime in a quantum gravity 
theory has been explored for many years, hrst in the context of black hole thermodynamics 
[30] and the Bekenstein bound [31] and later in more general holographic contexts, culmi¬ 
nating in the covariant entropy bound, he., the Bousso bound [32, 33]. 

The statement of the Bousso bound is the following: given a spacelike surface B of area 
A, draw the orthogonal null congruence on the surface and choose a direction in which the 
null generators have non-positive expansion. Let the null geodesics terminate at caustics. 
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singularities, or whenever the expansion becomes positive. The null hypersurface swept 
out by these null geodesics is called the lightsheet. Then the entropy S going through the 
lightsheet is less than A/AG. 

Let our spacelike surface B he a 2-ball of area ^4 on a spacelike slice of AdS and choose 
as the lightsheet the ingoing future-directed null congruence. This lightsheet will sweep out 
the entire interior of B and will terminate at a caustic at the center of B. Since the system 
is static, the entropy S passing through this lightsheet is the entropy of the system on B, 
which by the Bousso bound satishes 

sm < T . (22) 

It is natural to cast the Bousso bound as a constraint on the dimension of the bulk Hilbert 
space. As argued in Ref. [34], the thermodynamic entropy of a system about which we only 
know the boundary area A is just the logarithm of the dimension of the true Hilbert space of 
the bulk region in question (as opposed to the naive Hilbert space in quantum held theory), 
which the Bousso bound implies is less than A/AG} As such, if we denote the Hilbert space 
of B by l-isi let us replace Eq. (22) with the slightly more concrete statement 

In dim < —. (23) 

4Gr 

B. A MERA version of the Bousso Bound 

Our aim is to compute both sides of the inequality (23) using the MERA. For this 
calculation, it is instructive to change our parametrization of the hyperbolic plane from 
coordinates {x,z), which take values in the half-plane 2 ; > 0, to coordinates {p,0), which 
take values in a disk 0<p<l,0<6*< 27r. Embeddings of the MERA in a disk are often 
depicted in the literature, e.g., [38]; here we make this coordinate transformation explicit, 
since we wish to carefully study the geometric properties of the MERA. 

To begin, consider a MERA consisting of a single tree that contains a hnite number of 

layers m. This situation is illustrated in Fig. 6a for k = 2 and m = 4. Note that such a 

^ Moreover, it is known that there exists an asymptotically-AdS bulk configuration that saturates the Bousso 
bound, namely, the BTZ black hole [35, 36], which further implies that IndimRe in fact equals A/AG. 
However, we will not need this stronger assertion in what follows. A similar but unrelated result equating 
the area of a region with its entanglement entropy in vacuum was obtained in Ref. [37]. 
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MERA begins with a top-level tensor at the level that seeds the rest of the MERA in 
the IR. 





FIG. 6. (a) A k = 2 MERA consisting of m = 4 layers and with periodic boundary conditions, (b) 
the corresponding embedding in (x, z) coordinates, and (c) the embedding in (p, 6) coordinates. 


The base of the MERA is made np of sites. Withont loss of generality, let us locate the 
leftmost site of the base of the MERA at x = 0, so that the UV-most sites sit at coordinates 
(x, z) = {na, a), where n = 0,1, 2,..., (fc"* — 1) as shown in Fig. 6b. Let us also assume 
periodic boundary conditions for this MERA and hence identify x = 0 and x = 
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Next, define the coordinates {p,0) as follows: 


P 

e 


k^a 


271 , 


X 


k^a 


In these coordinates, the metric reads 


( 24 ) 


ds^ 


L 2 

(i^ 



(25) 


cf. Eq. (3). This embedding of the MERA is shown in Fig. 6c; the top-level tensor always 
sits at p = 0 and the lower layers of the MERA are equally spaced on circles of radii 
1/2, 3/4, 7/8, ... that are centered at p = 0. 


More generally, one could construct a top-level tensor that has T legs, each of which begets 
a tree of sites. In this case, a; = 0 and x = Tk"^~^a are identihed, so one should dehne the 
angular variable as 6 = 2'Kx/{Tk'^~^a). The metric (25) is correspondingly modihed and 
reads 


ds^ 


(l-p)2 + F V27rJ 


(26) 


This situation is depicted in Fig. 7. (If T = /c, however, then it is not necessary to introduce 
any new structure in addition to the disentanglers and isometries that were already discussed, 
he., one may take the top-level tensor to be one of the isometries.) 


We may immediately compute the right-hand side of Eq. (23). Let the ball B be centered 
about p = 0, and suppose B contains the top-level tensor, the sites at the top tensor’s legs, 
and then the hrst Ns generations of the MERA emanating from these sites, as indicated in 
Fig. 7. The boundary of is a circle at constant p, so its circumference according to the 
MERA is A = Tk^'^L. As such, we may write 


A _ Tk^BL _ Tk^Bc 
AG ~ AG ~ 6 


(27) 


where in the second equality we used the Brown-Henneaux relation, Eq. (10). 

How one evaluates the left-hand side of Eq. (23) using the MERA is not as immediate. 
Recall that 77^ is the Hilbert space of bulk states. The MERA, however, does not directly 
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FIG. 7. Disk parametrization of the Poincare patch of AdS in which a MERA has been embedded. 
The top tensor of the MERA shown has T = 6. The shaded region is a ball B, which is this case 
contains Ag = 1 generation. 


prescribe the quantum-gravitational state in the bulk; it is not by itself a bulk-boundary 
dictionary. As we mentioned in Sec. IIB, the minimal assumption that one can make is 
to posit the existence of a bulk Hilbert space factor Vbuik associated with each MERA site 
that is not located at the top tensor. To keep the assignment general, we assign a factor 
Vt to the top tensor. The dimensionality of each I4uik factor should be the same in order 
to be consistent with the symmetries of the hyperbolic plane. The assumption of a Hilbert 
space factor at every MERA site is minimal in the sense that it introduces no new structure 
into the MERA. A true AdS/MERA correspondence should dictate how states in the bulk 
Hilbert space are related to boundary states. However, for our analysis, it is enough to 
simply postulate the existence of the bulk Hilbert space factors Vbuik and Vt, each of which 
may be thought of as localized to an AdS-scale patch corresponding to the associated MERA 
site. 


In addition to the site at the top tensor, the number of regular MERA sites that the ball 
B contains is given by 


Nb 


/^JVb+1 _ 1 
^( k-l 


(28) 


i=0 ' ^ 

As such, the Hilbert space of bulk states restricted to B is T-Lb = (Vbuik)®'^® ® Vp. Next, 
suppose that dim Vbuik = X and that dimVr = Xt, where, like X) X and xt are some hxed, 
Ag-independent numbers. Then dim^fg = Xt(X'^®)- Note that one would expect x and x 
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to have a very specific relationship in a true bulk/boundary correspondence, the nature of 
which will be explored later in this section. Combining Eqs. (27) and (28), the dimensionality 
of T-Lb is upper bounded as follows: 


IndimT/B < tv:: 

4G 


T I —^^— 1 In X + In XT < 


- 


k-1 


6 


(29) 


After isolating c in Eq. (29) and using the result of Eq. (19), we find that 


In Y / — 1 1 

c = 12rj f(k)^ > 6 ( + 7 ^ IhXt ) • 


In k \k^'3{k — 1 ) 


Tk^B 


(30) 


Next, let us consider this inequality in the limit of large Ag. A motivation for this limit is 
the fact that the natural scale of validity of an AdS/MERA correspondence is super-AdS, as 
was established in Sec. III. Moreover, by virtue of its definition, there is always an ambiguity 
of order the AdS scale in the radius of the ball B. That is, the region in AdS denoted by 
B is only well-defined in the MERA if B is large compared to the AdS scale L. Taking the 
limit of large Ag, Eq. (30) reduces to 




kink 

2{k-l) 



(31) 


By using the bound on rjfXk) given by Eq. (18), we arrive at a constraint on k, x, and x^ 


k^lnk /lnx\ ^ 

2(k - 1 ):: ftox/ “ 


(32) 


In principle, the above inequality could be satisfied for any k, provided that the dimension 
X of the factors Vbuik can be arbitrarily chosen with respect to the bond dimension x- 
However, the essence of holography, that the bulk and boundary are dual descriptions of the 
same degrees of freedom and therefore have isomorphic Hilbert spaces [5] , implies a relation 
between x and x- Namely, for a MERA with a total of A levels of sites in the bulk strictly 
between the UV-most level and the top-level tensor, the number of bulk sites A/buik that are 
not located at the top tensor is given by Eq. (28) with Ag = A, and the number of sites in 
the boundary description is A/boundary = Tk^~^^. The bulk Hilbert space thus has dimension 
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and the boundary Hilbert space has dimension y^boundary_ Equating^ the dimension 
of the bulk and boundary Hilbert spaces then yields 


tax _ 1 Z'_ taXxA ^ j 

Iny A/bulk V Inx / 


(33) 


where we took the limit of N large, consistent with Eq. (31) and in keeping with the ex¬ 
pectation that the UV cutoff be parametrically close to the boundary at p = 1. Putting 
together Eqs. (32) and (33), we obtain a constraint on k alone: 


k'^lnk ^ 
2{k-l) - 


(34) 


This constraint cannot be satished for any allowed value of the rescaling factor k, which 
must be an integer greater than or equal to 2. We thus learn that a conventional MERA 
cannot yield a consistent AdS/MERA correspondence. The MERA cannot simultaneously 
reproduce AdS geodesics, respect the Ryu-Takayanagi relation, and (using the only con¬ 
struction for the bulk Hilbert space available to the MERA by itself) satisfy the Bousso 
bound. That is, there exists no choice of MERA parameters that can faithfully reproduce 
geometry, holographic properties, and bulk physics. 

If we relax this bound and, instead of Eq. (18), only observe the weaker, natural bounds 
7] < 1 and f{k)<k — las discussed at the end of Sec. IV C, the constraint (34) is corre¬ 
spondingly modihed: 


k\nk 

2(4-1) - 


(35) 


In contrast to Eq. (34), this latter bound can be satished, but only for k = 2,3, or 4. As such, 
other AdS/tensor network correspondences, in which the ancillae are perhaps entangled and 
therefore do not describe a conventional MERA, are not ruled out. Note that we never 
needed to compute bulk entanglement entropy explicitly — and therefore did not need to 
treat separately the possibility of entanglement among ancillae — because we cast the Bousso 
bound as a constraint on the size of the bulk Hilbert space itself. The appearance of p in 
Eq. (31) corresponds to entanglement in the boundary theory as computed by the tensor 
network; Eqs. (31) and (33) still apply. 

^ We recognize that there are other proposals [12, 39] that do not require an exact equivalence between 
the bulk and boundary Hilbert spaces, but, even in these cases, there is the requirement of an exact 
equivalence between the logical qubits on the boundary with the Hilbert space of the bulk. 
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VI. CONCLUSION 


The notion of emergence of spacetime based on a correspondence between AdS and a 
tensor network akin to AdS/CFT is a tantalizing one. A necessary step in snch a program 
is the evalnation and comparison of calculable quantities on both sides of the duality. In 
this work, we have subjected the proposed AdS/MERA correspondence to such scrutiny. To 
summarize, let us restate our three main findings: 

1. In matching the discrete graph geometry of the MERA to the continuous geometry of 
a spatial slice of AdS, we demonstrated that the MERA describes geometry only on 
scales larger than the AdS radius. Concretely, as shown in Sec. Ill, the proper length 
assigned to the spacing between adjacent sites in the MERA lattice must be the AdS 
scale. 

2. By requiring that the entropy of a set of boundary sites in the MERA — whose 
computation is a discrete realization of the Ryu-Takayanagi formula — be equal to 
the GET ground state entropy of the same boundary region in the thermodynamic 
limit, we obtained a constraint on the parameters that describe a MERA in terms of 
the GET central charge [Eqs. (20) and (21)], which implies that the bond dimension 
X must be exponentially large in the ratio of the AdS scale to the Planck scale. 

3. In the natural construction of a bulk Hilbert space ("Hbuik) using the MERA, we used 
the Bousso bound to constrain the dimension of 'Rbuik- When combined with our 
previous results, we found that any strict AdS/MERA correspondence cannot satisfy 
the resulting constraint, Eq. (34). Upon relaxing the definition of the MERA or 
allowing for additional structure, however, we obtained a looser constraint, Eq. (35), 
which may not rule out some other AdS/tensor network correspondences. 

In particular, more general correspondences between AdS and MERA-like tensor networks, 
in which we allow the ancillae to be entangled when reproducing the GFT ground state [and 
for which Eq. (35) applies in place of Eq. (34)] are not ruled out by our bounds, provided 
that the rescaling factor k = 2, 3, or 4. Further, it is interesting to note that our bounds 
extend to states other than the vacuum that are described by a MERA. One such example, 
namely, states at hnite temperature dual to black holes in AdS, is discussed in App. B below. 
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While the consistency conditions that we fonnd are specific to the MERA tensor network, 
many of the ideas and techniqnes that we nsed apply eqnally well to other tensor networks. 
In the EHM, for instance, the type of bnlk Hilbert space dimensionality argnments that we 
made based on the covariant entropy bonnd may be directly transferred to the EHM. The 
same stringent final constraints that we derived do not apply to the EHM, however, since it 
is nnclear to what extent the EHM reprodnces the Ryn-Takayanagi formnla (which renders 
the resnlts of Sec. IV inapplicable). Onr bnlk Hilbert space argnments similarly apply to 
the holographic error-correcting code proposal in Ref. [12], which fnrthermore pnrports to 
reprodnce a version of the Ryn-Takayanagi formnla. It is presently nnknown, however, 
whether the bonndary state of a holographic code can represent the gronnd state of a GET, 
so an identification of entropies similar to the identification Rmera = 5'cft, npon which onr 
bonndary entropy constraints so crncially depend, cannot yet be made. 

In closing, we have fonnd several consistency conditions that any AdS/MERA corre¬ 
spondence mnst satisfy. The totality of these constraints rnles ont the most straightforward 
constrnal of an AdS/MERA correspondence. Other interesting holographic correspondences 
that are described by tensor networks more general than the MERA and that respect all of 
onr bonnds may indeed be possible. Onr consistency conditions are nice validity checks for 
these correspondences when applicable and in other cases they may inspire similar consis¬ 
tency conditions. The program of identifying the emergence of spacetime from the bnilding 
blocks of qnantnm information is an ambitious one; stringent consistency conditions, such 
as those presented in this paper, are important for elucidating the subtleties in this quest 
and in providing guidance along the way. 
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Appendix A: Entropy bound for general MERAs 


Following the method presented in Ref. [22], let us compute an upper bound for the 
entanglement entropy of a region B consisting of £o sites in a MERA with rescaling factor 
k. We will use the notation of Ref. [22] throughout. 

First, recall the result from Ref. [22] that the entanglement entropy of a region consisting 
of Eq sites is bounded by 

5mera(4; B) < (O + N^^,) Iny. (Al) 

The quantity im' is the width of the causal cone at depth m! and N^, = Ylm=o is the 
total number of sites that are traced out along the boundary of the causal cone. In other 
words, N!^, is the number of bonds that are cut by the causal cone up to a depth m' {cf. 
Fig. 4). The quantity Iny is the maximum entanglement entropy that each site that is 
traced out could contribute to Smera{Eo;B). Note that Eq. (Al) holds for all m' > 0. 

The width of the causal cone for a given m' depends sensitively on the structure of the 
MERA. In particular, the number of sites that are traced out at each renormalization step 
depends on the choice of disentanglers, as well as how they are connected to the isometries. 
Eor instance, in a MERA with a rescaling factor k, any given disentangler could have any¬ 
where from 2 up to fc incoming and outgoing legs. (It should be reasonable to require that 
any disentangler can have no more than k incoming and k outgoing legs so that it straddles 
no more than two isometries.) It is thus clear that the number of bonds that one cuts when 
drawing a causal cone, and hence the entanglement entropy of the region subtended by that 
causal cone, depends on the choice of disentanglers and connectivity. 

Nevertheless, we can compute an upper bound for S'mera(^o;-B) by considering a worst- 
case scenario for the number of bonds cut by the causal cone. We begin by asking: What 
is the largest number of bonds that a causal cone could cut in one renormalization step 
at a depth m'l The layout of disentanglers and isometries that produces this situation 
is shown at one side of a causal cone in Fig. 8. If the causal cone at the bottom of the 
renormalization step incorporates a single bond that goes into a disentangler accepting k 
bonds, then the causal cone must cut the other k — 1 bonds entering the disentangler. Then 
if this disentangler is arranged so that its leftmost outgoing bond is the hrst bond to enter 
an isometry from the right, the causal cone must cut the other k — 1 bonds entering the 
isometry. If this arrangement is mirrored on the other side of the causal cone, we see that 


29 


4(fc — 1) bonds are cut by the causal cone in this renormalization step, he., = 4(fc — 1). 



FIG. 8. Left side of a causal cone that cuts the maximum possible number of bonds over the course 
of one renormalization step. The rectangles are disentanglers that accept k bonds as input and the 
triangles are isometries that coarse-grain k bonds into one. The causal cone is the shaded region. 
If this situation is mirrored on the right side of the causal cone, then 4(A; — 1) bonds are cut in this 
renormalization step. 


Recall that for any hnite io, after a hxed number of renormalization steps, the width of 
the causal cone remains constant for any further coarse-grainings. The depth at which this 
occurs is called the crossover scale and is denoted by m. Therefore, the causal cone will cut 
the largest possible number of bonds when the arrangement described above and depicted 
in Fig. 8 occurs at every step up until the crossover scale. Then, by Eq. (Al), the entropy 
bound is given by 

5'mera(4; B ) < {Ifa + 4(fc - l)m) In y, (A2) 

where is the width of the causal cone at the crossover scale. 

For any given causal cone in a MERA with scale factor fc > 2, the maximum number of 
additional sites the causal cone can pick up at some level vn! is 4(/c — 1). Therefore, for a 
causal cone that contains 1^' sites at depth m', the number of sites in the causal cone after 
one renormalization step ^m'+i < \{^m' + 4(/c — 1))/A:] < tm'Ik + 5. Applying the relation 
recursively, we hnd that the number of sites at any layer m' < m is bounded. 



m=l 


(A3) 


Setting m' = rh, it trivially follows that the crossover scale obeys rh < log;, io. Furthermore, 
we notice that this is the scale at which the entanglement entropy is minimized if we trace 
over the remaining sites. In other words, the number of bonds cut by going deeper into 
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the renormalization direction is no less than the bonds cut horizontally, so 4(fc — 1) > 
Applying the bounds for fh and Ifh on Eq. (A2), we arrive at an upper bound on S'mera(^o; B) 
for a fc-to-one MERA, 

5'mera(4;5) < 4(fc - 1)(1 + logfc4) Iny. (A4) 

When is parametrically large, we neglect the 0{1) contribution to the bound on 
*S'mera(4; B), which yields Eq. (13). 


Appendix B: BTZ Black Holes and Thermal States in AdS/MERA 


Thus far, we have found constraints on the structure of a MERA that can describe GET 
states dual to the AdSs vacuum. One might ask whether these results extend to other 
constructions that exist in three-dimensional gravity. Although pure gravity in AdSa has 
no local or propagating degrees of freedom, there exist interesting non-perturbative objects, 
namely, BTZ black holes [35]. In this appendix, we extend our constraints on boundary 
entanglement entropy to these objects. 

The non-rotating, uncharged BTZ black hole solution is given in Schwarzschild coordi¬ 
nates by 

(r2 _ 2\ t2 

ds' = + —i—dr^ + , (Bl) 

with a horizon at r = r_|_. Noting that Euclidean time is compactihed by identifying r ~ r-|- 
27rL^/r_|_, the horizon temperature of the black hole is given by T = r^j^'nL?. Additionally, 
the Bekenstein-Hawking entropy of the black hole is 


Area 7rr+ 


(B2) 


Let us now consider applying a MERA with rescaling factor k and bond dimension x 

to a GET at a hnite temperature, where instead of minimizing the energy of the boundary 

state, one minimizes the free energy. In the GET, turning on a temperature introduces a 

scale, going as the inverse temperature, which screens long-range correlations. Thus, the 

^ Alternatively, we can see this from a heuristic argument by noting that the crossover scale is the scale 
at which the causal cone has a constant width for further coarse-grainings, i.e., {im + 4(fc — l))/fc ~ fm- 
Therefore, fm ^ 4 < 4(/c — 1). 
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state will have classical correlations in addition to entanglement and the effect of a finite 
temperature on the entanglement entropy is the appearance of an extensive contribution. 
As one runs the MERA and coarse-grains, the thermal correlations that cannot be removed 
become more relevant. The MERA, which is unable to remove the extensive contribution, 
truncates at a level with multiple sites. The schematic entanglement renormalization process 
is illustrated in Fig. 9. The state at the top level effectively factorizes, where each factor 
appears maximally mixed [10, 21]. A tractable realization of this tensor network structure 
recently appeared in Ref. [40], which found a MERA representation of a thermal state. 


P P P P 



T 


^exp 



FIG. 9. The MERA, when applied to a thermal CFT state exp(—R cft/T), where Z = 
tr(exp(—R cft/T)), truncates after a finite number of layers. The boundary state at the top of the 
truncated MERA effectively factorizes into a product of maximally mixed states p = I jx- 


Keeping in mind that the holographic dual of a finite-temperature state in the CFT 
is a black hole in AdS, where the temperature of the CFT corresponds to the Hawking 
temperature of the black hole, we note that the truncated MERA is suggestive of a black 
hole horizon [10]. If the MERA is to be interpreted as a discretization of the geometry, then 
the geometry has ended at some scale. Also, as we approach the horizon, the amount of 
Hawking radiation that we see increases and the temperature measured by an observer at 
the horizon diverges. The density matrix of some system in the infinite-temperature limit 
is given by the product of a maximally mixed state at each site, just like the state at the 
top of the MERA. It is important to note that, as was pointed out in Ref. [40], in order 
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to reproduce the correct thermal spectrum of eigenvalues, a small amount of entanglement 


must be present between the sites at the horizon. If the bond dimension were taken to be 


inhnite, then the sites at the horizon truly would factorize. But for a hnite bond dimension, 
one should really think of the horizon as a high-temperature state, with sites effectively 
factorized. 

For small regions on the boundary, the length of the subtending bulk geodesic is subexten- 
sive and so the Ryu-Takayanagi formula maintains that the boundary region’s entanglement 
entropy is subextensive as well. However, if we consider a large enough region on the bound¬ 
ary, the geodesic will begin to probe the horizon of the black hole. The geodesic will run 
along the black hole horizon and pick up an extensive contribution to the entropy. We 
consider a boundary theory living on a lattice consisting of rib sites, with total system coor¬ 
dinate length Xgys = n\ya. In the limit as r approaches the boundary in the metric (Bl), we 
see that Txsys = ‘r+/L, as was pointed out in Refs. [6, 29]. We further note that this implies 
that the system coordinate size is of order AdS radius, x^ys = 27rL. 

Let us now view the MERA of Fig. 9 as a discretization of a BTZ spacetime and repeat 
the analysis of Sec. III. In this discretization, the layers of the MERA lie along circles of 
hxed radius r in the coordinates of Eq. (Bl). Again, we ask what proper length Li separates 
sites in any given layer of the MERA. 

First, note that a path at hxed tq that subtends an angle 0o has proper length ro0o- At 
the boundary of the MERA, we consider a region dehned by 0 < 0 < 0o = 27ra;o/xsys, where 
Xq is the coordinate length of the interval, consisting of Iq lattice sites. The boundary of the 
MERA is at a hxed radius r = rb. Naturally, the boundary radius rb can be interpreted as 
a UV cutoh and is related to the lattice spacing a by rb = L^/a [6]. By equating the proper 
distance of the region in the MERA, f'o-hi, with that at the boundary of the BTZ spacetime, 
we hnd the proper length between horizontal bonds to be Li = L. 

With the foresight that the top of the MERA is suggestive of a black hole horizon with 
proper length 27rr_|_, the number of sites at the hnal layer is therefore rib = 27rr_|_/L. This 
further tells us that the MERA truncates after a hnite number of layers m, given by 



(B3) 


This coincides with the conclusion in Refs. [40, 41] that the MERA representation of a 
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thermal state is obtained after C>(log^(l/T)) iterations of coarse-graining. 

Now consider a region B on the boundary consisting of sites and for which the cor¬ 
responding geodesic contains a segment running along the BTZ horizon. The subextensive 
contribution to the entropy in the MERA is exactly as before, in which we pick up at most 
Iny from each bond we cut with the causal cone of the region B. Furthermore, we will 
now pick up an extensive contribution from the horizon, where the number of horizon sites 
within the causal cone is and each such site in the product state on the horizon contributes 
maximally to the entropy by an amount Iny. Combining the contributions, we hnd 

5'mera(5) = 4r7B/s(fc)logfc lnx-h41nx. (B4) 

Recall that the entanglement entropy of a single interval B of coordinate length xq in a 
CFT at hnite temperature [27] is given, up to a non-universal constant, by 

S'cft(5) = 

where xq is much smaller than the total system size Xgys- The standard held-theoretic 
derivation of the above entropy is done by computing the Euclidean path integral on an 
n-sheeted Riemann surface and analytically continuing to hnd the von Neumann entropy. 
The same result can be derived by computing geodesic lengths on spatial slices of BTZ 
spacetimes and making use of the Ryu-Takayanagi formula. 

When T —0 in Eq. (B5), we recover the usual result (9). In the T —>■ cx) limit, the 
von Neumann entropy gives the usual thermal entropy as entanglement vanishes. Taking 
Txq S> 1, the leading and subleading contributions to the entanglement entropy are 

5'cft = -vtxqT H — In-— , (B6) 

3 3 27raT ’ ^ ’ 

where the hrst term is the thermal entropy for the region B. 

Now let us consider a hnite-temperature CFT that is dual to a BTZ black hole with 
horizon temperature T = r_|_/27rL^. In terms of geometric MERA parameters, we hnd that 
Eq. (B6) becomes 

^CFT = ^4 + ^mlnfc. (B7) 

D 3 
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Here we used the fact that t'h = xqt^/L'^ as well as Eq. (B3), where we note that m can 
also be written as log^(£b/f'h)- The result (B7) coincides precisely with the extensive and 
subextensive contributions calculated using the MERA in Eq. (B4) provided that c/ In y; ~ 
0{1). Therefore, we hnd that the truncated MERA correctly captures the entanglement 
structure of thermal GET states and their dual BTZ spacetimes. These conclusions are in 
agreement with those in Refs. [20, 41]. 

As a check of the claim that c and Iny should be of the same order, we can compare 
the horizon entropy given by the contribution from the sites at the hnal layer with the 
Bekenstein-Hawking entropy (B2) of a BTZ black hole. There are Uh sites comprising the 
horizon, each with Hilbert space dimension y. The system is in the inhnite-temperature limit 
— and hence described by a maximally mixed density matrix, with entropy contribution 
In y from each site — so 

5'horizon = Uh In y . (B8) 

Making use of the Brown-Henneaux relation and requiring that the entropy (B8) coincide 
with the Beckenstein-Hawking entropy, we again hnd that c/lny ~ 0{1). More specihcally, 
taking the counting to be precise, we hnd that 


c/ Iny = 6, 


(B9) 


which is qualitatively in agreement with the previous conclusion (20) that the Hilbert space 
dimension must be exponentially large in c. 


With this relation, the extensive terms in Eqs. (B4) and (B7) agree precisely. Further 
identifying the subextensive terms, we hnd = (ln/c)/2. If we then impose the 

constraint (18), we hnd that 


kink 

2(4-1) - ^ ' 


(BIO) 


This last inequality exactly reproduces Eq. (35) and thus constrains k to be 2, 3, or 4. 
Interestingly, we have found the weaker of the two bounds derived in Sec. V, without needing 
to consider the Bousso bound. 


As desired, the truncated MERA computation of entanglement entropy agrees with the 
expected entanglement entropy given by the application of the Ryu-Takayanagi formula to 
the length of the minimal surface in a BTZ spacetime. The fact that the results of matching 
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boundary entanglement entropy given in Sec. IV further hold in BTZ spacetimes might not 
be too surprising given that such spacetimes are quotients of pure AdSs. 
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